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ABSTRACT 

We present the first detection of orbital motion in tlie cataclysmic variable 
EX Hydrae based on X-ray data from the Chandra X-ray Observatory. The 
large collecting area of the telescope and the high resolution of the HETG 
spectrometers allow for an unprecedented velocity accuracy of ~15 km s^^ in 
the X-ray wavelength regime. We find an emission line velocity amplitude of 
58.2±3.7 km s~^ and infer a white dwarf mass of 0.49 ± 0.13 Mq, in good agree- 
ment with previous studies using optical, ultraviolet, and far ultraviolet data. 

Subject headings: novae, cataclysmic variables — stars: individual (EX Hydrae) — 
techniques: spectroscopic — x-rays: stars 



1. Introduction 

Magnetic cataclysmic variables (CVs) are semi-detached, interacting binaries consisting 
of a late-type star that overflows its Roche-lobe and a white dwarf with a magnetic field suffi- 
ciently strong to influence the accretion of matter originating from the late-type companion. 
The magnetic field directs the accretion fiow away from the orbital plane and onto both of the 
white dwarf's magnetic poles. Accretion behaves differently in the two classes of magnetic 
cataclysmic variables. In polars the stars are tidally locked; this "static" orientation results 
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in an accreting flow from the Lagrangian point LI of the binary along the magnetic field 
lines of the white dwarf onto its magnetic pole(s). In intermediate polars (IPs), where the 
white dwarf is spinning faster than the orbital period, matter also accretes onto the white 
dwarf along the magnetic field lines; however, the matter first accumulates in an accretion 
disk surrounding the white dwarf. In both polars and IPs the accreting matter approaches 
the poles of the white dwarf supersonically and so passes through a stand-off shock before 
coohng and setthng onto the white dwarf (Aizu 1973; Kylafis & Lamb 1982). 

The cataclysmic variable EX Hydrae belongs to the IP class. It consists of a white 
dwarf with a spin period of 67 minutes and has an orbital period of 98 minutes. Both the 
spin period and binary period can be identified in its light curves, from optical (e.g., Vogt, 
Krzeminski, & Sterken 1980; Sterken et al. 1983), UV (e.g., Mauche 1999; Belle et al. 2003), 
and X-ray emission (e.g., Watson, Sherrington, & Jameson 1978; Allan, Hellier, & Beardmore 
1998). Furthermore, the inchnation of the EX Hya system is such that its hght curve shows 
a partial eclipse of the lower accretion pole of the white dwarf by the companion (Allan et 
al. 1998). 

Optical and UV spectroscopy of the EX Hya system has detected the orbital motion 
associated with the white dwarf (e.g., Breysacher & Vogt 1980; HeUier et al. 1987; Mauche 
1999; Belle et al. 2003). This motion, combined with the orbital velocity of the companion 
and the presence of the eclipse, allows for a mass determination of the white dwarf inde- 
pendent of models of the accretion column. In this letter we use Chandra High-Energy 
Transmission Grating (HETG) data for EX Hya and show that its orbital motion can be 
detected with high resolution X-ray spectra. 

2. Observation and Reduction 

EX Hya was observed by Chandra using the HETG in combination with the Advanced 
CCD Imaging Spectrometer in its spectroscopy layout (ACIS-S) on 2000 May 18 for 60 ks. 
The observation was continuous and covers ~10 orbital revolutions of the binary system and 
~15 white dwarf revolutions. 

We reduced the data using the Chandra Interactive Analysis of Observations (CIAO) 
software package^ making only two departures from the standard reduction: 1) we turned 
randomization off (i.e., we set rand_pix_size = 0.0 in the tg_resolve_events tool) to min- 
imize any artificial line broadening, and 2) we applied a solar system barycentric correction 
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so that the event times are in Bary centric Dynamical Time instead of spacecraft time. The 
reduced HETG EX Hya spectrum, as shown in Figure 1 (see also Mauche 2002; Mukai et 
al. 2003), reveals emission lines from H- and He-like ions, as well as from Fe L-shell ions, 
superimposed on a weak continuum. Fe XVII and Fe XXII density diagnostics reported by 
Mauche, Liedahl, & Fournier (2001, 2003) indicate densities of order Ug ~ 10^^ cm~^. 

3. Velocity signature 

To search for velocity signatures in the EX Hya HETG data we used the following 
procedure: 1) identify the spectral lines in the first-order High Energy Grating (HEG) and 
Medium Energy Grating (MEG) spectra; 2) divide the observation into time segments; 3) 
add the spectral lines in each time segment to form a composite line profile (CLP) for each 
grating arm in order to boost the signal-to-noise ratio; 4) measure the line centroid of the 
CLP; and 5) look for variations of the line centroid with orbital or white dwarf spin phase. 

3.1. Line Identification and Wavelength Calibration 

Table 1 lists the isolated spectral lines, i.e., lines that do not have another strong line 
within two times the instrumental FWHM, with a signal-to-noise ratio larger than 10. We 
obtained the reference wavelengths for these lines from the Astrophysical Plasma Emission 
Database (APED; Smith et al. 2001). Reference wavelengths with uncertainty estimates 
from APED are taken from quantum electrodynamic theory for H-like ions (Ericsson 1977) 
and He-hke ions (Drake 1988) and from laboratory measurements for Fe L-shell fines (Brown 
et al. 1998, 2002). Wavelengtfi uncertainties are for individual lines only and do not include 
tfie effects of dielectronic satellite fines that are associated with many of these strong lines. 
Such blending satellite lines contribute less than a few percent to the strength of the emission 
feature. 

We defined unblended lines as eitfier single lines or multiplets tfiat are unresolved and 
from tfie same ion. For unresolved lines from tfie same ion we used tfie emissivity-weigfited 
wavelengtfis. Tfie He-like Mg and Si forbidden and intercombination lines are unresolved by 
MEG, and since tfieir tfieoretical line flux ratios depend on electron density and pfiotoexciting 
radiation, we cannot weigfit tfieir wavelengtfis by tfie low-density emissivities tabulated in 
APED. Tfius we excluded tfieses fines from tfie CLP. 

To confidently add spectral lines it is important to verify tfie wavelengtfi calibration and 
to correct for any systematic offsets between tfie observed line centroids and tfie reference 
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wavelengths from the APED database. To fit the emission hnes we first determined the 
contribution of the continuum emission to the spectrum. A global fit to line-free regions 
of HEG and MEG spectra using the Astrophysical Plasma Emission Code model (APEC, 
version 1.3, Smith ct al. 2001) allows a self-consistent treatment of the continuum for all 
lines. We defined the line-free regions as any part of the spectrum greater that 0.1 A in 
width that is devoid of spectral lines with emissivities larger than 10~^^ ph cm^ s~^ (this is 
~1/1000 of brightest line in the spectrum) in the APED database. We found a good fit for a 
thermal continuum (consisting of bremsstrahlung, radiative recombination, and two-photon 
emission) for a temperature of 16±1 keV (see Figure 1), in good agreement with the shock 
temperature of 15.4l2;6 keV reported by Fujimoto Sz Ishida (1997) for an Aizu (1973) model 
fitted to the ASCA spectrum of EX Hya, and the maximum temperature of 20 keV reported 
by Mukai et al. (2003) for a cooling fiow model fitted to the same Chandra HETG data. A 
two-temperature continuum model did not significantly lower the of the fit. This is not 
surprising since the highest temperatures dominate the emission measure distribution in an 
accretion column cooling mainly by means of thermal bremsstrahlung. 

To determine if any offsets between the observed wavelengths and their reference values 
exist we performed the following fit. Wc set the distances between each consecutive pair 
of emission lines equal to their separation based on the APED database and only fit the 
normalization of each fine and the overall wavelength offset. Each fine is modeled using a 
Gaussian with a FWHM equal to the instrumental FWHM (0.012 A for HEG and 0.023 A 
for MEG) . This resulted in the following offsets for each of the orders (predicted — observed 
wavelength): -0.0012(2), -0.0020(2), -0.0008(2), and 0.0004(2) A for HEG-1, HEG+1, 
MEG— 1, and MEG-l-1, respectively. Although the offsets are small, they do constitute ~10% 
of the HEG FWHM and ~5% of the MEG FWHM, or about 30-50 km s'^ (at 10 A), and 
need to be corrected for. 

Even more important than the absolute and relative wavelength calibration is its sta- 
bility during an observation; this appears to be extremly stable ( Chandra Proposers Obser- 
vatory Guide^). Furthermore, the rotation and orbital periods of the EX Hya system are 
not relevant timescales for the satellite (dither period is 707.1 sec in the pitch direction and 
1000 sec in the yaw direction). 



^http: / / cxc.harvard.edu/ proposer/POG/index.html 
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3.2. Binciry 



To look for changing velocities associated with the binary period in the Chandra HETG 
data for EX Hya we first phase folded the spectra (using the ephemeris of Hellier & Sproats 
1992) and divided the observation into phase intervals centered at (pgg = 0.05, 0.10, 0.95, 1.00 
with a width of A^gg = 0.25, where ^gg indicates the phase in the binary orbit. Note that 
the phase intervals are overlapping. 

Unfortunately, by splitting the observation into smaller phase intervals none of the 
spectral fines have sufficient signal-to-noise ratio to confidently measure changes in their 
wavelengths. We therefore added all fines in each interval to accumulate enough counts to 
accurately measure the central wavelength. We accomplished this by the following procedure: 
1) For each grating first-order arm we determined the count- weighted average wavelength of 
the lines in Table 1, A = CjAicf,*/ where q is the number of counts in the spectrum 

for a segment 2xFWHM wide centered on Aicf.i (from Table 1) minus the number of counts 
predicted by the model continuum (see §3.1). 2) Each line in the spectrum was shifted by an 
amount equal to its reference wavelength, Arcf.i and by the wavelength calibration correction 
(see §3.1). 3) The wavelength-binned shifted spectrum was then transformed to a velocity- 
binned spectrum according to vun,! = c(Abin,i — Aref,i)/Aref,i. 4) The shifted velocity-binned 
profile for each spectral line was added to form one line profile. However, since the bin size 
for each line was different we chose to rebin each t'bin,? to the bin size corresponding to A: 
^^bin = c(Abin — A)/A. 5) Finally, we measured the central velocity for this composite line 
profile (CLP) . 

Figure 2 shows the CLPs for binary phases ^gg = 0.00, 0.25, 0.50, and 0.75. The change 
in CLP centroid is clearly visible showing how powerful the combination of a large effective 
area and high spectral resolution can be in the X-ray regime. 

The radial velocity shifts as a function of ^gg for the combined orders MEG— 1 plus 
MEG-l-1 are shown in the top panel of Figure 3. We present the MEG data rather than the 
HEG data because they have a factor of 2-10 times more counts per spectral fine, which 
improves our centroiding accuracy even with the lower resolution of the MEG instrument 
(see also §5). Results from the HEG and the individual MEG orders are consistent but with 
lower significance. The figure clearly shows that there is a systematic motion of the CLP 
with a period equal to that of the binary. The solid line represents the following fit to the 
data points: 




(1) 



-6- 



where 7 = 1.3±2.3 km s~\ K = 58.2±3.7 km is the velocity amplitude, 0o = 0.48±0.01, 

ui is fixed to 1.0, and A = 0.25. Since the data points overlap, and are thus correlated, we 
calculated fo^^ each of the five sets of independent data points (indicated in Figure 3). 
The psr degree of freedom (dof) for the independent data sets range from x^/(dof = 1) 
= 0.1 to 1.0. It is important to note that the value of x^ for a small number of dof can 
be different from the value one expects to indicate a good fit in the case of a large number 
of dof, i.e., x^/dof = 1.0. The important statistic is Q(x^| dof), which is the probability 
that the observed chi-squared will exceed the value x^. Acceptable values for Q, i.e., the 
3(7 confidence interval, range from ~0.0015 to ~0.9985. For the above values of xV^of, Q 
ranges from 0.3 to 0.8. 

We also fitted the data with 1) a constant velocity and 2) a sine curve with the period 
not forced to the binary phase, i.e., u not fixed to 1.0 in Eq. 1. The former fit yields 
7 = -2.8 ± 2.3 km s~\ xV(dof = 3) ranges from 13 to 19, and Q f« 7.0 x 10-^\ while the 
latter fit yields a period of = 1.01 ± 0.04 times the binary period, and values of 7, K, 
and 00 are indistinguishable from the fit described above. The latter fit has zero degrees of 
freedom and hence a meaningless x^- 

The above velocity amplitude and phase offset, K — 58.2 ± 3.7 km s~^, and 0o = 
0.48 ± 0.01 are in good agreement with those reported by Brcysacher & Vogt (1980, K = 

68 ± 9 km s"\ 0o = 0.5), Cowley, Hutchings, & Crampton (1981, fsT = 90 ± 28 km 

00 = 0.55), Gilliland (1982, K = 58 ± 9.3 km s"\ 0o = 0.52), Hellier et al. (1987, K = 

69 ±9 km s-\ 0o = 0.53 ±0.03), Mauchc (1999, X = 85 ±9 km s-\ 0o = 0.54 ±0.02, using 
FUV lines), and Belle et al. (2003, K = 59.6 ±2.6 km 0o = 0.48 ±0.05, using UV lines). 

3.3. White Dwarf Spin 

The bottom panel of Figure 3 shows the CLP radial velocity versus the white dwarf spin 
period using bins of size one-third the spin period (i.e., A = 0.333). The velocities have been 
corrected for the average binary motion for each data point based on the fit described in §3.2. 
To provide upper limits on the rotation velocity of the white dwarf and the infall velocities 
in the accretion column, we fitted sine curves to these data. We assumed that the upper 
accretion column is pointing away from the observer at white dwarf phase (consistent with 
both the accretion curtain model [Rosen et al. 1988] and the accretion pole occultation model 
[Allan et al. 1998]) and that the cooler part of the lower accretion column, where most of the 
spectral lines are formed, is not visible throughout most of the spin cycle (Allan et al. 1998). 
The dotted curve in the figure represents the best-fit rotation velocity of 9.2 ± 5.0 km s^^ 
(00 = 0.50) and the dashed curve represents the best-fit infall velocity of 7.9 ± 5.4 km s~^ 



- 7- 



(00 = 0.25). The 3(7 upper limits on the rotation and infall velocities are thus ~25 km s ^. 

4. White Dwarf Mass 

Combining the orbital motion of the white dwarf obtained in §3.2 with the orbital motion 
of the companion (Vande Putte et al. 2003, K2 — 360 ± 35 km s~^) and the echpse duration 
(Mukai et al. 1998) allows us to determine the white dwarf mass Mi, companion mass M2, 
and inchnation i of the EX Hya system. 

Here we followed the approach taken by Beuermann et al. (2003) using a computer 
program kindly provided to us by K. Beuermann. Assuming that the emission originates 
at or near the pole of the white dwarf, the iterative method yields Mi = 0.49 ± 0.13 M©, 
M2 = 0.078 ± 0.014 M0, i = 77?2 ± 0?6. The white dwarf mass Mi corresponds to a white 
dwarf radius Ri = 1.0 ± 0.2 x 10^ cm (Wood 1995). If the emission originates above the 
surface of the white dwarf the inclination of the system decreases ~ 0?1 for every 0.1 x Ri 
in elevation above the surface. These results agree within the errors with those found by 
Beuermann et al. (2003) and hence do not change any of the conclusions reached in that 
paper. 

Finally, we note that our estimate of the white dwarf mass, which is based on the radial 
velocity of the emission lines in the Chandra HETG spectrum, agrees with that of Fujimoto 
& Ishida (1997) (Mi = 0.48lo:o6 M©), which is based on the H- to He-hke emission hue 
intensity ratios observed in the ASCA spectrum. Our "classical" radial velocity result is 
applicable only to high-inclination magnetic CV binaries, whereas theirs is applicable to a 
wide range of polars and IPs, although it relies on theoretical models of accretion column. 

5. Discussion 

Our radial velocity measurement in EX Hya pushes the instrumental limit of the Chan- 
dra grating instruments. The MEG instrument has an effective area large enough to accu- 
mulate between 3000 and 4000 counts in the CLP so that we can rehably find the centroid 
to an accuracy of ~0.5 mA, corresponding to an accuracy in velocity space of 10-15 km s"^ 
(at ~10 A). This is consistent with the statistical uncertainty for measuring Gaussian line 
centroids derived by Landman, Roussel-Dupre, & Tanigawa (1982): 

5, ^a ^-j {AxfJ'A, (2) 
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where 5"^ is the statistical uncertainty in /i (the hne center) , a is the instrumental resolution 
(FWHM/2.3548), (Aa;)o is the bin size expressed in terms of a, and A is the average error 
expressed in terms of the peak amplitude. For our MEG CLP this equation yields 0.5 mA, 
as observed, where we have taken A = 2/\/P, with P the number of counts in the peak bin 
of the CLP. 

Unfortunately the effective area of the HEG instrument is significantly smaller (a factor 
of 3 at 10 A) and does not allow us to improve over the MEG despite the better instrumental 
resolution (S*^ = 0.5 mAfor HEG, similar to MEG). This means that for this EX Hya data 
set we obtain similar velocity accuracies for HEG and MEG. To fully exploit the HEG 
resolution, i.e., gain a factor of 2 in velocity accuracy, a significantly longer observation 
would be required. 

More velocity resolution is desirable to detect velocity changes induced by the rotation 
of the white dwarf around it spin axis. Such velocity signatures would provide information 
on the location and velocity of the emitting plasma and on the location of the magnetic 
pole. The lack of a pronounced velocity variation at the spin period of the white dwarf in 
this HETG observation is expected. To see this, note that the free-fall velocity onto the 
white dwarf vg — (2G'Mi/i?i)^/^ fa 3600 km s~^, and the immediate post-shock velocity 

V — ffr/4 900 km s~^. Also note that the majority of the spectral lines in the HETG 
spectrum are formed at temperatures of less than 2 kcV, i.e., 0.1 times the post-shock 
temperature T = SGMiiriYifi/SkRi ^ 20 kcV (/x = 0.615 is the mean molecular weight of 
a completely ionized plasma of solar composition). From Figures 1 and 2 of Aizu (1973) 
we estimate the velocity of the 2 keV gas to be 0.1 w ~ 90 km s^^. Assuming that the 
accretion columns arc perpendicular to the orbital plane, the projected post-shock velocity 

V cosi ~ 18 km s~^. Furthermore, the maximum rotational velocity of the white dwarf 
i^Tot — 27rit!i/67min ^ 15 km s~^. Both velocities are smaller than the upper limits derived 
in §3.3. 



6. Conclusions 

We have confirmed the orbital velocity of the white dwarf in the cataclysmic variable 
EX Hya using the Chandra HETG. This measurement is the first detection of orbital motion 
in a CV binary system based on X-ray spectroscopy. By adding various spectral lines visible 
in the EX Hya spectrum enough signal-to-noise ratio can be accumulated to obtain velocities 
accurate to 10-15 km s^^. No new conclusions are drawn on the white dwarf mass (0.49±0.13 
Mo). 
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Fig. 1. — Spectrum for the combined first orders of the MEG instrument. Smooth curve 
represents the one-temperature thermal continuum model fit to the spectrum. 
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Fig. 2. — Composite line profiles (CLPs) based on the first-order MEG data for the non- 
blended lines of EX Hya for binary phases 0.00, 0.25, 0.50, and 0.75 (shifted in the vertical 
direction by 200, 400, and 720 counts, respectively, to separate the profiles). The average 
error on the CLP bins is depicted in the last bin of each CLP. Dotted line profiles represent the 
instrumental profile of the CLP (normalized at the peak) at a constant velocity of 1.3 km s~^. 
Number of counts in the 0.75 profile is lower because this binary phase corresponds to the 
bulge absorption feature in the EX Hya light curve. 
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Fig. 3. — Top panel: CLP radial velocity versus the binary period based on the first-order 
MEG data. Sohd curve represents the sinusoidal fit to the data. The five sets of independent 
points have been indicated by different symbols (see §3.2). Bottom panel: CLP radial velocity 
(corrected for the binary motion) versus the white dwarf spin period. Dotted curve represents 
the best-fit velocity curve corresponding to the rotation velocity of the white dwarf and the 
dashed line represents the best-fit velocity curve corresponding to the infall velocity in the 
accretion column (see §3.3). 
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Table 1. Strong EX Hydrae HETG spectral lines 



Ion 


Upper level 


Lower level 


Wavelength'' 

rAi 


Fe XXVI 




Is 25 , 


1.780 


sXVI 


2r> In q/9 

1/^,0/^ 


Is 251/2 


4.729 


sXV 


ls2p iPi 


ls2 ISo 


5.039 


Si XIV 


2n /o -I/O 


Is ^Si/a 


6.182 


Si XIII 


ls2p 


ls2 ifi-,) 


6.648 


Si XIIP 


ls2p 2 


ls2 'So 


6.687 


Fe XXIV 


4p ^^^1/2,3/2 


2s 251/2 


7.989 


MgXII 


2p 2^1/2,3/2 


Is 25i/2 


8.421 


MgXI 


ls2p iPi 


ls2 'So 


9.169 


MgXIb 


ls2p 3Pi,2 


ls2 i5o 


9.231 


NeX 




Is 25i/2 


10.239 


Fc XXIV 


■id 2/35/2 


2s 25i/2 


11.176 


Fe XXIV<= 


3s 25j/2 


2p 2P3/2 


11.432 


Fe XXIII 


2sZd 'D2 


2s2p 'Pi 


11.736 


Ne X^ 


2p ^-Pi/2,3/2 


Is 2Si/2 


12.132 


Ne IX 


ls2p iPi 


ls2 'So 


13.447 


Fe XVIII 


2pi/22p^/23d6/2 


2p5 2^3/2 


14.208 




2p4(iD)3d 2D5/2 


2p'' 2p3/2 




Fe XVII 


2p5(2p)3d ipi 


2p6 'So 


15.014 


Fc XVII 


2p"(2p)3(i ^^Di 


2p« i5o 


15.261 


Fc XVIIP 


2p4(3p)3s 2P3/2 


2p5 "P3/2 


16.005 


Fe XVII 


2p5(2p)3s iPi 


2^6 i5o 


16.780 


Fe XVII 


2p5(2p)3s 3pj 


2p6 i5o 


17.051 


OVIII 


2p ^^1/2, 3/2 


Is "Si/2 


18.968 



"AH wavelength errors arc better than 0.004 A. 

'^The unresolved He-like forbidden and intercombination lines 
for Si and Mg are not considered in the velocity analysis. 

<=Blend consists of the Fe XXIV, Fe XVIII 2p4(3p)4d 2^5/2 
to 2p5 2p3^2, and Fe XXII 2s2pi/23p3/2 to 2s22p lines. 
We adopted the wavelength for Fe XXIV since most flux seems 
to be contained in this line. 

"^Blend consists of the Ne X Ly a doublet, Fe XVII 
2p5(2p)4d ipj to 2p6 i5(,, and Fe XXII 2p{^P)3d "D^/^ to 
2p2 '^D^i2 lines. 

''Blend consists of the Fe XVIII and O VIII 3p 2Pi/2 3/2 to 
Is 25i/2 lines. 



